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1.  Introduction 


The  effect  of  bottom  topography  on  the  flow  of  rotating  and/or 
stratified  fluids  has  recleved  a great  deal  of  attention  in  the  recent 
past.  The  experiments  by  Sir  G.I.  Taylor  (1921)  in  which  the  so-called 
Taylor  column  was  first  observed  in  the  laboratory  were  the  genesis  of 
these  efforts.  McCartney  (1975)  presents  a table  delineating  the  most 
Important  theoretical  studies  on  Taylor-columns  including  the  parameter 
restrictions  for  the  various  analyses. 

Theoretical  studies  for  the  homogeneous  incompressible  case  include 
those  of  Jacobs  (1964),  Ingersoll  (1969)  and  Vaziri  and  Boyer  (1971). 

The  first  two  of  these  investigations  were  restricted  to  inviscid  flows: 
i.e.  E^<<Ro  where  E and  Ro  are  the  Ekman  and  Rossby  numbers  respectively. 
Unfortunately  laboratory  experiments  on  this  class  of  flows  are  generally 
restricted  to  parameter  ranges  in  which  viscous  effects  are  important; 
i.e.,  E^'X/Ro.  The  Vaziri  and  Boyer  study  developed  a viscous  model  in 
which  E^'V/Ro  and  in  addition  presented  a series  of  laboratory  experiments 
which  were  in  good  agreement  with  the  theory  advanced.  The  experiments 
utilized  a rotating  water  tunnel  (analogy  with  a wind  tunnel  in  a non- 
rotating  frame)  which  was  found  to  provide  improved  observational  cap- 
abilities to  the  towing  experiments  first  conducted  by  Taylor  and  later 
considered  by  Hide  and  Ibbetson  (1966) . 

More  recently  Davies  (1972)  has  extended  the  towing  experiments  to 
include  the  effects  of  stratification.  Whereas,  the  homogeneous  incom- 
pressible experiments  were  generally  restricted  by  the  Taylor- Proudmen 
i- 

I 

I theopem  which'  implies  a horizontal  velocity  field  independent  of  the 

I 

i distance  along  the  axis  of  rotation,  these  stratified  experiments  allowed 
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changes  in  the  velocity  field  with  depth.  The  Davies'  experiments  again 
were  in  the  general  range  Ro'^'E^  so  that  viscous  effects  were  presumably 
Important;  specifically  Davies  utilized  the  following  set  of  parameters: 
Ro  - 8.8(10)”^  and  E = 4.0(10)"^. 

Hogg  (1973)  in  a theoretical  study  followed  up  the  Davies'  experi- 
ments with  an  investigation  of  the  effect  of  stratification  on  Taylor 
columns.  This  study  neglected  viscosity  and  investigated  three  ranges 
of  stratification:  S«Ro;  S'''Ro;  and  SMDCl)  where  S is  a stratification 
parameter. 

McCartney  (1975),  neglecting  viscous  effects,  presented  an  analysis 
of  a two-layered  rotating  flow  on  a beta-plane;  he  also  conducted  a 
series  of  towing  experiments  for  the  single-layer  case.  In  another 
beta-plane  study  Janowitz (1975)discussed  an  inviscid  solution  for  flow 
of  a continuously  stratified  fluid  over  shallow  topographies. 

In  the  present  study  the  flow  of  a continuously  stratified  rotating 
fluid  over  shallow  topographies  is  considered.  The  initial  objective  of 
the  study  is  to  develop  a theoretical  model  for  a rotating  stratified 
flow  which  can  be  examined  in  the  laboratory  utilizing  the  rotating  water 
tunnel  concept  advanced  earlier  by  one  of  us;  e.g.  see  Boyer-Ridge  (1971), 
Boyer-Step  (1971),  or  Vaziri  and  Boyer  (1971).  That  is  an  experiment  is 
to  be  developed  in  which  a stably  stratified  fluid  is  to  be  in  uniform 
motion  relative  to  a rotating  observer.  Topographic  features  are  then  to 
be  placed  in  the  moving  stream  with  the  objective  of  determining  the 
resulting  flow  field  and  comparing  same  with  the  theoretical  studies. 

. In  the  hbmogeneous  incompressible  case  the  laboratory  experiment 
consisted  of  the  flow  through  a channel  of  rectangular  cross-section. 
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The  flow  outside  of  the  boundary  layers  on  the  channel  walls  could  be 
shown  to  be  uniform  and  geostrophlc.  Because  there  were  no  free  surfaces. 
It  could  be  shown  that  the  centrifugal  terms  In  the  equations  of  motion 
were  unimportant  and  thus,  that  the  location  of  the  channel  with  respect 
to  the  rotation  axis  was  not  Important. 

When  stratification  is  Introduced,  however,  the  position  of  the 
rotation  axis  with  regard  to  the  tunnel  test-section  will  be  an  Important 
consideration.  In  the  analysis  given  in  Sections  2 and  3 we  assume  an 
Idealized  physical  system  in  which  the  flow  is  forced  to  move  azimuthally 
about  the  axis  of  rotation  in  an  annulus  whose  upper  and  lower  bounding 
surfaces  are  parabolic  (Figure  1).  These  parabolas  are  the  equilibrium 
surfaces  of  a given  stably-stratified  fluid  system  at  rest  with  respect 
to  the  rotating  system.  A specified  topographic  feature  is  located  on 
the  lower  surface  and  it  is  desired  to  determine  the  resulting  flow  field 
for  a specified  upstream  shear  flow. 

In  section  2 the  governing  equations  for  the  resulting  flow  are 
derived  for  a range  of  parameters  considered  attainable  in  the  laboratory. 
In  section  3 specific  solutions  for  the  flow  over  a long  ridge  of  constant 
cross-section  and  an  axisymmetric  "bump"  are  given. 

In  Section  4 progress  on  the  laboratory  experiments  is  given. 


2.  Governing  Equations 


Consider  the  physical  system  sketched  In  Figure  1 a stably  strati- 
fied fluid  Is  confined  In  an  annulus  defined  by  two  parabolic  surfaces 
separated  by  a distance  H and  vertical  walls  at  r = R — j-  and  r = R + — . 
The  annulus  Is  rotating  at  a constant  angular  velocity,  about  a verti- 
cal axis  and  the  fluid,  away  from  solid  boundaries.  Is  assumed  to  be  In 
uniform  motion  with  speed  dependent  on  the  distance  from  the  lower  sur- 
face In  the  azimuthal  direction.  As  discussed  earlier  the  Ekman  layers 
on  the  parabolic  bounding  surfaces  combined  with  the  constraints  of  the 
vertical  walls  would  tend  to  alter  the  stratification.  It  Is  assumed 
here,  however,  that  the  Ekman  layers  can  be  controlled  l.e.,  by  removing 
fluid  from  the  corner  regions  so  that  a stable  stratification  can  be  main- 
tained. 

Define  an  Intrinsic  coordinate  system  in  which  x is  in  the  azimuthal 
direction,  y is  along  the  surface  with  posistive  being  toward  the  axis  of 
rotation  and  z is  normal  to  the  surface.  A topographic  feature  h(x,  y) 
of  characteristic  horizontal  length,  L,  and  vertical  scale,  h,  is  located 
in  the  vicinity  of  the  origin.  It  is  desired  to  determine  the  motion  in 
the  vicinity  of  the  topography. 

Define  the  density  field  as 

P = Pq  + Pgfz)  + P(x,  y,  z) 

where  p^  is  the  specified  density  at  z = 0 assuming  no  motion,  P^Cz)  is 
the  specified  vertical  distribution  assuming  no  motion  and  p is  the 
density  field  compatible  with  the  motion  field.  Define  Ap^  = Pg(H)  - p^  (0) 
as  the  characteristic  magnitude  of  p^. 
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Also  define  the  pressure  field  as 


p = Pq(z)  + p(x,  y,  z) 


where 


- (Pq  + Pg(z))  G k 


2 A 2 

Here  G = (g  + u R and  p is  the  pressure  associated  with  no  motion. 

1*  U)^l< 

Taking  H to  be  of  the  order  of  L and  restricting  to  •5<<1  and  <<1, 

K 8 

the  steady-state  equations  of  momentum,  diffusion  and  conservation  of 
mass  can  be  written  as 

(p  + p + P)  -^  = -Vp  - 2(p  + p + p)  (o)  X v)-  p(G  + O(a)^L))  k 

0 s Dt  os 

*y  0 

+ p 0(u)  L)j  + pV  V 


dPg  2-.  ‘^^Ps 

(v.V)p  + w = kV^p  + K — -V" 

dz  dz2 


(2.1) 


7.pv  = (v.7)p  + w -jj—  + (p^  + p^  + p)  7.V  = 0 

2 

where  v(u,  v,  w)  is  the  velocity;  V = -I—  i + j + ^ + 0(-^  and  ; 

dv  dv  d R D 


— = V.7;  p the  viscosity  and  k the  mass  diffusivity. 
Define  the  following  dimensionless  quantities: 

V.  P . _ E 


X*  =^;V*=-X;  3*  = 


2o)Up  G 
0 


-1 


; p 


. Q * = 

2a;Up  L’  ^s  fip 


(2.2) 


where  we  anticipate  the  density  and  pressure  field  to  be  determined  by  a 
balance  of  gravitational  and  coriolis  effects. 

Introducing  (2.2)  into  (2.1)  one  obtains: 


D /I  ^^3  \DVj.t72Dv  ^ /I.  ^^3  2a)U^ 

Rod  + “ ^3^  dT  dF  = ■ ~ — )(toxv) 

0 0 

' 2 2 
- ^ (1  + 0(^^i))k+  pO(^)j  + E7^F 
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— « If  2^  —1  2» 

Ro(v.V)p  + B 77  ” = in:  B 7 + Ea  7 p 


d2 


O IT  ^P  1 n 

2a)U  /■-  . s svi+vij.  s , 2uU  _ - » 

(v.V)p  + -—  — Po  + —F'  p)  7.V  = 0 

O p d2  p s G 
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(2.3) 


2 2 

where  Ro  = U/2a)L  is  the  Rossby  number;  Fr  = U /GL  the  square  of  the 


2 2 

Froude  number;  E = v/2toL  the  Ekmdn  number;  B 


^'PsG,,  2 


P L 
o 


/4o)  the  stratifi- 


cation parameter;  and  a = v/ic  the  Prandtl  number. 

Relations  (2.3)  are  five  equations  in  the  unknowns(u,v,(i),p,p)  and 
will  be  solved  in  approximate  form  subject  to  the  no  slip  condition  on 
the  upper  and  lower  bounding  surfaces.  It  will  be  assumed  that  the 
lateral  bounding  surfaces  will  not  effect  the  motion  in  the  vicinity  of 
the  topography.  Upstream  of  the  topographic  feature  it  is  assumed  that 
the  motion  is  given  by  U(z).  Downstream  we  require  the  velocity  to  be 
azimuthal. 

The  reduction  of  (2.3)  to  a tractable  form  is  guided  by  the  objective 
of  deriving  solutions  in  a region  of  parameter  space  that  can  be  examined 
in  the  laboratory.  In  this  regard  we  consider  typical  parameters  used 
in  the  similar  homogeneous  experiments  conducted  by  Vaziri  and  Boyer 
(1971) and  further  use  stratification  values  that  could  be  attainable  from 
a saline  stratified  water  experiment;  e.g. 


U = 0.1  cm/sec 
(1)  = 0.1  rad/sec 


G Iv  g = 980  cm/sec 


V = 0.01  cm^/sec 

K = 1.6(10)  ^ cm^/sec 
3 


p^  = 1.01  gm/cm 


L = 2.5  cm 
H = 2.5  cm 
h = 0.25  cm 
R = 100  cm 


Ap.  = 0.01  'gm/cm 
s 


■s 


ii 
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from  whence 


E - 8.0(10) 
Ro  = 0.2 
= 97.0 

Fr^  = 4.1(10) 


-X 

0 


1.6(10) 


H 

L 


1.0 


h 

L 


0.1 


I = 2.5(10)"^ 


We  now  assume  the  following  relationships  for  the  dimensionless 
parameters: 


E«1 


Ro  = kE^(k  unity) 

r.2 

B unity 

Fr^  'v-  (a  > 0) 


a ^ 'Vi  E^  (y  > 0) 


H 


— unity 


^(x,y)^  h^(x,y)  E^  (h^'v-  unity)| 


>(2.4) 


It  is  assumed  that  v,  p and  p can  be  expanded  in  power  series  in 


i 

E^:  i.e. 


- - o i 

v=VE  +VE^+... 
o 1 

p = S^E°  + S^E^  + . . . \ 

p = P E°  + P,  E^  ' ■ ■ 

^ o 1 


(2.5) 


The  zeroth  order  is  dictated  by  the  method  of  non-dimensionalization  while 
the  first  order  is  suggested  by  the  character  of  the  Ekman  layers  on  the 
parabolic  surfaces.  Substituting  (2. 5)  into  (2.3)  and  utilizing  restrictions 
(2.4)  one  obtains: 


First  Order 
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0 - - P + V 
ox  o 

0 - - P ^ U 
oy  o 

0 - - p _ s 

02  . O 


0 


(2.6) 


U + V + W =0 
ox  oy  02  I 

Second  Order 

'v 

- k(U  U + V U ) = - P,  + V, 
o ox  o oy  lx  1 

k(U  V + V V ) = - P,  - U, 
o ox  o oy  ly  1 

0 - - P,  - 1 R Z£_  U ~ S,  \ 

4 L Ro  ° ^ / (2.7) 

dp„ 

k(U  S + V S ) + b2  W = 0 
o ox  o oy  dr  1 

U,  + V,  + W,  =0 
lx  ly  l2 

^ J 

Cross-differentiating  the  first  two  of  (2.7)  and  utiliring  the  last 


of  (2.7)  gives  an  expression  for  W,  in  terms  of  zeroth  order  terms. 

This  can  then  be  equated  to  W,  as  determined  from  the  fourth  of  (2.7). 

Iz 

Utilizing  the  third  of  (2.6)  then  one  determines  the  following  relation 


for  P 


(2.8) 


If  we  assume  a linear  profile  for  p„  (i.e.,  Pq  = - „ intro- 

duce  a new  vertical  coordinate  ^ define  by  C = B(— )^z,  (2.8)  reduces  to 

n 


J(P_^,  72  P^)  . 0 


(2.9) 


where 


= -L  + ^ + -5- 

^3  ^2.22 

3x  3y  3z 


Under  the  above  scaling,  the  boundary  layers  on  the  parabolic  sur- 
faces are,  to  lowest  order, Ekman  layers.  The  boundary  conditions  on 


these  surfaces  are  the  same  as  those  given  in  Vaziri  and  Boyer  (1971). 
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Along  the  lower  surface  the  Ekman  layer  solutions  specify  the  vertical 
velocity  at  the  outer  edge  of  the  layer  (i.e.  z » 0 of  the  interior)  as 


*1  “vTr 


1 ? + J(P  , h ) 


(2.10) 


while  along  the  upper  surface  (z  = — ) 

W « — ih 

1 z Z o 


(2.11) 


2 2 

2 d a 

where  ^2  “ — 2 — 2‘  vertical  components  must  match  the  interior 

3x  3y 

values  specified  by  the  fourth  of  (2.7)  from  whence  one  obtains 


JT  ‘'oz>  ■ ^ ^ ■ 0 


(2.12) 


<rf  i "oz>  ■ » « = ■ ! 


(2.13) 


respectively.  The  boundary  conditions  far  upstream  and  downstream  are 


P *~1,  X and  P * 0,  x ■+■  + 

oy  ’ ox  ' 


(2.14) 


respectively.  At  large  lateral  distances  from  the  topographic  feature 


P = -1,  y -»-  + 
oy  - 


(2.15) 


The  task  is  thus,  to  solve  (2.9)  subject  to  (2.12)  through  (2.15). 
From  (2.9)  we  note  that 


7,P  = f(P  ) 

JO  o 


(2.16) 


Then  from  (2.14)  one  notes  that  V^P  far  upstream  is  zero  for  all  P and 

JO  o 

hence  f(P  ) is  identically  zero.  The  analytical  problem  thus  reduces  to 
o 

a solution  of  Laplace's  equation  subject  to  the  non-linear  set  of  bound- 
ary conditions  (2.12)  to  (2.15).  Let  us  now  consider  some  specific 


examples. 


3 


. Examples 
A.  Long  Ridge 

We  first  consider  a ridge  with  infinite  extent  in  the  y-direction 
and  having  a uniform  triangular  cross-section;  specifically  take 


h (x) 
o 


- f 


X < - i 
< X < 0 

0 < X < J 

X > J 


(3.1) 


Assuming  the  velocity  components  to  be  Independent  of  y,  the  formulation 
of  Section  2,  now  expressed  in  terms  of  the  velocity  components,  leads  to 
the  following  boundary  value  problem: 

Governing  Equations 


V + V =0 

oxx  L ozz 

o 


"lx,  * "izz  - 0 

D 


(3.2) 


Boundary  Conditions 


r 


Wj^(x,o) 


Wi(x,  l)  - 

I 

Wj^(x  - oo,  z) 


W^(x 


=>,  z) 


7T  '^ox 

X < -i 

-i  < X 1 0 

-o  V 

v2  ox 

0 < X 1 ♦ 

^ V 
'/2  ox 

X ^ i 

. V 

ox 

- 0 

V (x  z)  • 0 

O 

0 

V (x  -+•  »,  z)  * 0 

O 

(3.3) 


A 
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where  o = — E The  numerical  solution  of  (3.2)  subject  to  (3.3)  is 
straightforward  and  leads  to  and  velocity  profiles  as  exemplified 
by  Figures  2 and  3 respectively. 

For  flow  of  a homogeneous  fluid  over  a similar  ridge,  Boyer-Ridge 

(1971)  found  V to  be  negative  for  the  entire  flow  field.  It  is  noted 
o 

here  that  is  positive  in  the  vicinity  of  the  ridge  on  its  upstream 
side.  Similar  results  were  found  by  Janowitz  (1975)  in  his  inviscid 
study.  One  should  also  note  the  decreased  horizontal  deflection  of  the 
flow  field  at  large  vertical  distances  from  the  toporgaphic  feature. 

This  is  in  constrast  to  the  homogeneous  case  for  which  these  deflections 
are  indipendent  of  depth. 

B*  Axisymmetric  Bump 

Next  consider  an  axisymmetric  bump  given  by  the  relation 

where  h,  is  the  bump  height,  R its  radius,  and  r the  radial  coordinate  in 
d 

the  x-y  plane.  Laplaces  equation  in  was  thus  solved  numerically  sub- 
ject to  the  boundary  conditions  (2.12)  to  (2.15).  Figures  4A,  4B,  and 
4C  give  the  horizontal  streamline  pattern  for  the  flow  at  the  "bottom", 
"middle"  and  "upper"  surfaces  of  the  flow  field.  We  note  the  rapid  de- 
crease of  the  effect  of  the  topographic  feature  as  the  vertical  distance 
from  the  feature  increases.  We  also  note  that  the  flow  field  at  the 
lower  levels  is  reminiscent  of  the  flow  in  the  homogeneous  case  consid- 
ered by  Boyer  and  Vaziri. 

Figures  5A,  5B,  and  5C  give  the  vertical  velocity  field  at  the 

"bottom",  "middle"  and  "upper"  surfaces  respectively.  The  "low"  pressure 
regions  on  the  diagram  depict  areas  in  which  the  vertical  velocity  com- 
ponent is  negative  (i.e.  downward)  while  the  high  pressure  regions  depict 


Figure  2.  Ridge  flow  V profiles;  E = 0.00078,  Ro  = 0.039 
b2  = 0.97,  H?L  = 0.75  and  h/L  = 0.063. 


Figure  3.  Ridge  flow  profiles;  same  parameters  as 
Figure  2. 
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a positive  vertical  velocity  component.  The  profile  along  the  lower 
surfaces  Is  substantially  different  then  the  two-cell  pattern  found  by 
Vazlrl  and  Boyer  for  the  homogeneous  case.  Vazlrl  and  Boyer  found  a 
negative  cell  that  generally  covered  the  topography;  l.e.,  a suction  of 
fluid  toward  the  topographic  feature.  Upward  flow  was  found  in  a cell 
downstream  and  to  the  right  of  the  topography.  The  principal  new  aspect 
for  the  stratified  case  is  the  strong  positive  cell,  upstream  and  to  the 
right  of  the  topography.  Figures  5B  and  5C  depict  the  rapid  decrease 
in  strength  of  the  vertical  velocity  component  with  elevation. 


rrrr 


Figure  4C  Horizontal  streamline  pattern;  same  parameters  as  4A 
except  z = 0.50. 


Figure  5A  Vertical  velocity  field;  same  parameters  as  4A. 


A.  Laboratory  Experiments 
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A series  of  laboratory  experiments  was  initiated  to  establish  a 
stratified  flow  in  a rotating  system  which  could  be  utilized  to  test  the 
theory  advanced  above.  To  date  a system  has  been  developed  which  is 
capable  of  providing  a uniform  flow  of  a two-layer  fluid  system  in  a 
rotating  frame. 

Figure  6 is  a schematic  diagram  of  the  system.  The  basic  system 
consists  of  simultaneously  moving  a lighter  fluid  (Tank  A)  and  a heavier  fluid 
(Tank  B-salt  stratified)  from  a non-rotating  frame  through  a test  section 
in  the  rotating  frame  and  hence  to  a disposal  tank  in  the  non-rotating 
frame.  Fluid  A,  for  example  starts  in  Supply  Tank  A and  thence  is  piped 
by  a pump  through  a block  valve,  a throttling  valve  and  a flow  meter, 
into  a reservoir  mounted  on  the  rotating  table.  The  fluid  then  flows 
under  gravity  through  a flow  meter,  a block  valve  and  a throttling  valve 
into  the  test  section  of  the  rotating  tunnel.  The  fluid  is  then  removed 
from  the  tunnel  to  a disposal  tank  located  an  the  rotating  table. 

The  above  system  is  now  undergoing  a final  series  of  tests  to 
establish  a uniform  flow  in  the  test  section.  The  system  has  worked  well 
in  that  the  fluids  move  through  the  tunnel  witha distinct  interface. 
Refinement  in  the  horizontal  excursions  of  the  flow  from  that  of  a uniform 
rectilinear  free-stream  are  required  before  tests  with  particular  obstecles 
can  be  conducted. 
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Figure  6.  Stratified  Rotating  Water  Tunnel 
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